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In the early Universe, or near a supernova core, neutrino flavor evolution may be affected by 
coherent neutrino-neutrino scattering. We develop a microscopic picture of this phenomenon. We 
show that coherent scattering does not lead to the formation of entangled states in the neutrino 
ensemble and therefore the evolution of the system can always be described by a set of one-particle 
equations. We also show that the previously accepted formalism overcounts the neutrino interac- 
tion energy; the correct one-particle evolution equations for both active-active and active-sterile 
oscillations contain additional terms. These additional terms modify the index of refraction of the 
neutrino medium, but have no effect on oscillation physics. 
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I. INTRODUCTION 

It is well known that interaction with a medium mod- 
ifies neutrino dispersion relations or, equivalently, gives 
the medium a nontrivial refraction index for neutrinos. 
Because this effect is generically flavor-dependent, it can 
have a profound impact on neutrino flavor evolution. In- 
side the Sun and the Earth, the refraction effect arises as 
a result of neutrino interactions with electrons and nucle- 
ons. On the other hand, in the early Universe or near a 
supernova core, where the number density of neutrinos is 
sufficiently high, the refractive properties of the neutrino 
background itself are important 0, S S S 13 ■ This neu- 
trino "self-refraction" is the subject of the present work. 

Early studies @, 0, II H H of the effect treated the 
neutrino background analogously to the case of ordinary 
matter (electrons and nucleons), namely, assuming that 

• for each neutrino, one can write a single-particle 
evolution equation and take into account the effect 
of all other particles (including other neutrinos) 
by adding appropriate terms to the one-particle 
Hamiltonian; 

• the contribution of the background neutrinos to 
this Hamiltonian is diagonal in the flavor basis. 

With these assumptions, the evolution of each neutrino 
in the system was described by 



' dt 



« 



(1) 



where i and j label neutrinos in the system, H vac and 
i? ma t are the usual vacuum and ordinary matter Hamil- 
tonian terms and the last term is the sum over the contri- 
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butions of all background neutrinos, taken to be flavor- 
diagonal. 

These assumptions were critically reexamined by 
James Pantaleone H who made several important ob- 
servations regarding their validity. First, he reasoned 
that flavor evolution of neutrinos in a neutrino back- 
ground is in general a many-body phenomenon and it 
is a priori not obvious that the first assumption is jus- 
tified in all cases. Second, he demonstrated that, for a 
system of several active neutrinos flavors, even when the 
first assumption is justified, the second one is incompat- 
ible with the symmetry of the problem. Indeed, the low 
energy neutral current Hamiltonian |30j |. 



NO 




(2) 



where the currents j£ = VaH^v a and a is a flavor index, 
a = 1,..N, possesses a U(N) flavor symmetry, which 
must be respected by any effective description of the sys- 
tem. This requirement is not satisfied by the diagonal 
form for Hj/P used in the earlier studies. 

Pantaleone proposed a modified form for H v t , which 
contained non-zero off-diagonal terms, 
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In Eq. the wave function of the background neutrino 



is normalized such that J dV(\. 
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and the coefficient of proportionality A equals \/2Gf(1 — 
cos P^) 7 P 1 -- 1 ^ being the angle between the two neutrino 
momenta. For antincutrinos in the background, the form 
of the Hamiltonian is exactly the same, with the only 
difference that A has the opposite sign. As can be easily 
checked, Eq. J3J is indeed consistent with the U(2) flavor 
symmetry of the two-flavor system. 

The result © was also later obtained by Sigl and Raf- 
felt [HI and by McKellar and Thomson [n) in the frame- 
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work of more general analyses of the flavor evolution of a 
neutrino system, which took into account both refraction 
effects and collisions. Eqs. llll.'tl) have been used as a start- 
ing point for extensive studies of the neutrino evolution in 
the early Universe ^ElMJOM 13 El El El Hj and 
in supernovae 0, l22l l2.il I24L l25l | . The genera l prope rties 
of the solutions have been investigated in [51 El El E3 • 

Certain theoretical aspects of these equations, how- 
ever, have not received adequate attention in the litera- 
ture to date. One such aspect is the validity of a descrip- 
tion of the problem by a set of single-particle equations. 
The absence of quantum correlations (entanglement) at 
all times during the evolution of the system is a priori 
not obvious, and has not been proven. Indeed, in the 
derivations 0] and E3 ^ is stated as an assumption. In 
0, it is suggested that for quantum correlations not to 
form the system must obey certain physical conditions, 
for example, it should contain an incoherent mixture of 
mass eigenstates. This raises the question how general 
the results in Eqs. lll.'il) are and what physical criteria 
determine their breakdown. 

The second important aspect is the connection be- 
tween the elementary neutrino-neutrino scattering pro- 
cesses and the macroscopic description of refraction given 
in Eqs. (HIMIi . For refraction in ordinary matter, such 
a connection is very well established and provides the 
most straightforward derivation of the effect as an in- 
terference of many elementary scattering amplitudes. A 
similar treatment for the case of the neutrino background 
has not been given. 

In this paper we present a description of the flavor 
evolution of neutrinos in a neutrino background in terms 
of the elementary neutrino-neutrino scattering processes. 
This description provides a transparent, physical picture 
of the effect, and, at the same time, allows us to prove 
that quantum correlations in the neutrino ensemble are 
negligible and therefore the description in terms of one- 
particle equations is valid. These equations are obtained 
directly from our formalism and are compared to the ac- 
cepted results, Eqs. (11 IMP . 

We work in the regime in which (i) neutrino-neutrino 
interactions are described by the low-energy four-fermion 
Fermi coupling, (ii) the neutrino gas is non-degenerate, 
and (iii) incoherent scattering of neutrinos with other 
neutrinos and other particle species is negligible. The 
first condition implies that the neutrino center-of-mass 
energies are well below the weak scale. The second one 
is satisfied if a gas of neutrinos has a number density 
n v <C E 3 , where E is a typical neutrino energy. Finally, 
incoherent scattering is significant if the column density 
of the medium exceeds the inverse of the scattering cross 
section: d = J n(x)dx <; 1/cr, as happens, e.g., in the 
early Universe at temperatures larger than few MeV (see, 
e.g., 0). 

The refraction effects of the neutrino background are 
negligible if the coupling between a neutrino and the neu- 
trino gas is significantly smaller than the vacuum oscilla- 
tion Hamiltonian or than the coupling to ordinary mat- 
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FIG. 1: Neutrino beam traversing a slab of FCNC interacting 
matter. 



tcr. This condition thus depends on the neutrino energy, 
on the density of the neutrino background, and on the 
density of the ordinary matter. In the core of the Sun, 
where the number density of neutrinos is n v ~ 10 6 cm" 3 , 
neutrino-neutrino interaction is negligible for all relevant 
neutrino energies and oscillation parameters. In contrast, 
near a supernova core the neutrino density can be as high 
as n v ~ 10 31 cm" 3 and is comparable with the electron 
density there. Therefore, neutrino-induced refraction ef- 
fects are important and must be taken into account. 

The paper is organized as follows. In Sect.[n]we review 
the microscopical picture of the neutrino flavor evolution 
in normal matter and point out that a naive extension of 
this picture to the neutrino self-refraction leads to seem- 
ingly paradoxical results. In Sect. Hill we show how this 
picture should be constructed consistently by identify- 
ing states in the neutrino ensemble that amplify coher- 
ently. We also show that coherent scattering does not 
form entangled states. In Sect. II Vl we compare the one- 
particle evolution equations we obtain for active-active 
and active-sterile oscillations in the neutrino and antincu- 
trino backgrounds with the accepted results. In Sect. IVI 
we show how entanglement is effectively destroyed by the 
refraction phenomenon. Finally, Sect. lVll summarizes our 
conclusions. 



II. FORMULATION OF THE PROBLEM 
A. Conventional neutrino refraction: FCNC case 

We begin by briefly reviewing the physics of the re- 
fraction effect in "normal" matter. While there are many 
ways to derive the refraction properties in this case, fun- 
damentally, the relevant effect is the coherent interference 
of many elementary scattering events. 

Let us consider the problem of neutrino oscillations in 
a medium which possesses flavor changing neutral cur- 
rent (FCNC) interactions. Such interactions give rise to 
the nonzero off-diagonal terms in the neutrino evolution 
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Hamiltonian, 



H 



V2G 



Fn 2 



FCNC — 



const 



(4) 



where Gf is the Fermi constant and n 2 is the number 
density of scatterers in the medium. 

As a toy example, consider a beam of electron neutri- 
nos incident on a thin slab of matter of thickness L made 
of FCNC interacting particles, as illustrated in Fig. ^ 
Assume that the neutrino masses are sufficiently small so 
that the effects of vacuum oscillation can be neglected. 
The flavor conversion rate in the slab can then be found 
using the following straightforward physical argument. 
Let / be the amplitude for an electron neutrino to scat- 
ter as a muon neutrino in a given direction on a particle in 
the target. If the scattering amplitudes for different tar- 
get particles add up incoherently, the flux of muon neutri- 
nos in that direction is oc N s \f\ , where N s is the number 
of scatterers. In the case of forward scattering, however, 
the scattering amplitudes add up coherently and, hence, 
the forward flux of muon neutrinos is oc N 2 \f\ 2 . Indeed, 
in the small L limit Eq. gives 



P. 



FCNC 



e 2 (G F n 2 L) 2 /2 



(5) 



which has the form Py^ v c oc A^l/I 2 , since e oc /. No- 
tice that by choosing a small L limit we were able to 
ignore the secondary conversion effects in the slab, i.e., 
to assume that for all elementary scattering events the 
incident neutrinos are in the v e state. 

To summarize, for small enough L, the flavor conver- 
sion rate due to coherent FC scattering in the forward 
direction is proportional to the square of the modulus of 
the product of the elementary scattering amplitude and 
number of scatterers. This quadratic dependence on N a 
is what makes the coherent forward scattering important 
even when the incoherent scattering can be neglected. 

Notice that exactly the same arguments apply if one 
considers the usual flavor-diagonal matter term due to 
the electron background in a rotated basis, for instance, 
in the basis of vacuum mass eigenstates. In this basis, 
the matter Hamiltonian has off-diagonal terms, resulting 
in transitions between the vacuum mass eigenstates. 



B. Neutrino background: physical introduction 

We seek the same description for the case of neutrino 
background. Let us therefore modify the setup in Fig. ^ 
and replace the slab by a second neutrino beam, such 
that the neutrino momenta in the two beams are orthog- 
onal (see Fig.[5Jl. To keep the parallel between this case 
and the FCNC case, we will continue to refer to the orig- 
inal beam as "the beam" and to the second beam as "the 
background" . The neutrinos in each beam can be taken 
to be approximately monoenergetic |3ll| . We again as- 
sume that the neutrino masses are sufficiently small so 
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FIG. 2: Toy problem to illustrate neutrino flavor conversion 
in the neutrino background. 
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FIG. 3: Elementary scattering event that causes a change of 
the flavor composition of the beam 



that, although flavor superposition states could be cre- 
ated outside the intersection region, the effects of vacuum 
oscillation inside the intersection region can be neglected. 
Any flavor conversion that takes place in the system is 
therefore due to neutrino-neutrino interactions in the in- 
tersection region. 

Let us first compute the amount of flavor conversion 
in the beam using Eqs. The conversion is expected 

because of the presence of the off-diagonal terms in these 
equations. The result depends on the flavor composition 
of the background. If the background neutrinos are all 
in the same flavor state 



cos av e 



sin av H 



(0) 



and their density is n 2 , the Hamiltonian for the evolution 
of a beam neutrino takes the form 



H = 



V2G 



Fn 2 



const 



cos 2a 
sin 2a 



sin 2a 
- cos 2a , 



(7) 



4 



After traversing the intersection region, a neutrino in 
the beam will be converted to the state with the prob- 
ability 

P Ue ^ ~ sin 2 2a(G F n 2 L) 2 /2, (8) 

assuming as before that the size of the region L is small. 

The above descriptions of neutrino flavor conversion in 
the neutrino background and in the FCNC background 
are very similar in form. We will next show, however, 
that despite very similar appearances of the equations, 
the underlying physics in the two cases is different. 

First, we need to establish which elementary processes 
give rise to neutrino flavor conversion in the neutrino 
background. In the absence of vacuum oscillations, the 
only interactions in the problem are neutral current inter- 



actions between pairs of neutrinos. These are described 
by the Hamiltonian J5J), which conserves the total flavor 
of the system. The flavor composition of the beam there- 
fore changes only if some of the background neutrinos of 
different flavor scatter into the beam, i.e., if a neutrino 
from the background and a neutrino from the beam ex- 
change momenta, i/(k) + i/(p) — > v(p) + is(k). As observed 
in [t| , such events can add up coherently and give rise to 
the flavor off-diagonal entries in the oscillation Hamilto- 
nian. 

Let us consider such an elementary event, as depicted 
in Fig. Following we restrict scattering angles to 
the directions in which coherent interference occurs and 
write the Hamiltonian for an interacting pair of neutrinos 
in the form 



dt 



/ \Ve(k)Pe(p)) \ 
Wn{k)v e {p)) 

V K(fcK(p)> / 



V2Gi 
V 



(1 - cos/3) 



2 
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2 



/ MfcK(£)> \ 

\Vii{k)v e {p)) 



(9) 



where V is the normalization volume. For a state \S), 
which at t = is \v e fa), this equation formally has a 
solution 

= - K"e) + e - i5 Et W v v) + hvO ^ 

where 5E = 2v / 2Gf(1 — cos/3)/V. In practice, however, 
one only needs a small t expansion of Eq. 11011 : 

S(t) ~ \v e ^} -i5Et/2(\v e v^ + 1*^)) 

= (1 + ia)\v e v li } + ia\ufj,iy e }, (11) 

where a = — \/2Gf(1 — cos (3)t/V. Indeed, since the in- 
teraction time and the normalization volume are deter- 
mined by the size of the neutrino wave packet I, t ~ I and 
V ^ I 3 , the absolute value of a in any realistic physical 
situation is always much less than 1. This simply means 
that the interaction between two neutrinos is described 
by the lowest order Feynman diagram and second order 
scattering effects may be safely ignored. 

We now return to the situation depicted in Fig. [21 Wc 
are interested in the muon neutrino content of the beam 
after it crosses the interaction region. Using the Hamil- 
tonian in Eq. to describe an elementary scattering 
event u e + v x , we find 

S(t) ~ \v e v x ) +ia(\v x v e ) + \v e v x )), 
= (1 + ia)\v e Vx) + iacosa\v e v e ) 
+ ia sin a\u ' e ) . (12) 

Thus, for an elementary event, the amplitude of mea- 
suring the neutrino in the beam as is proportional 



to sin a. If one multiplies this amplitude by the num- 
ber of scatterers and squares, as was previously done 
in the FCNC example, one finds P Vti ^ v oc iV 2 sin 2 a, 
which is in clear conflict with the prediction of Eq. JSJ, 
P Ue ->v u oc W s 2 sin 2 2a. 



Which of the two results is right? At first sight, the 
first possibility may appear more plausible. Consider, 
for instance, what happens when an electron neutrino 
propagates through a background of muon neutrinos 
(a = 7r/2). Since muon neutrinos appear in the beam as 
a result of elementary exchanges between the beam and 
the background, one may think that for a pure muon neu- 
trino background the conversion rate in the beam should 
be maximal. The first result indeed has this behavior, 
while the second result predicts no conversion. One may 
therefore be tempted to conclude that the second result 
fails to describe the system in this limit and perhaps that 
this failure signals the general breakdown of the validity 
of the single-particle description. 



As will be shown in the next section, however, despite 
its seemingly paradoxical behavior, the second result is 
actually correct. The explanation of the paradox lies in 
the procedure of adding up elementary amplitudes, which 
in the case of the neutrino background is more subtle 
than in the case of ordinary matter. 
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III. MICROSCOPIC ANALYSIS OF THE 
EFFECT 

The key observation is that a neutrino-neutrino scat- 
tering event changes not only the state of the beam, but 
also the state of the background. This implies that the 
refraction by the neutrino background is intrinsically dif- 
ferent from the case of ordinary matter, and requires a 
specific analysis to establish what states in the neutrino 
ensemble are amplified coherently. 

We begin by recalling that an elementary scatter- 
ing event can only be amplified if the particles scat- 
ter "forward" , meaning that either particle momenta do 
not change (t-channel), or that they are exchanged (u- 
channcl) as in Fig. |21 The interaction is described by 
Eq. Ullfl . For the sake of clarity, in this section we con- 
sider the exchange diagrams only and omit the effects 
of the non-exchange diagrams. The latter produce an 
overall (flavor- independent) phase shift of the neutrino 
states, and do not affect the neutrino flavor conversion. 

We first consider a single electron neutrino in the beam 
interacting with several neutrinos in the background. 
Assume that these background neutrinos are all in the 
same flavor state v x , Eq. JfJJ), and that the state of the 
background initially is a product of single-particle states, 
X X X ■ ■ ■ X X ). The total system therefore is initially in the 
state \e)\xxx...xx) and, as a result of the interaction, 
evolves over a time step St according to 

\e)\xxx...xx) — ► \F) = \e)\xxx...xx) +ia\F\), (13) 

where \F{) is the state with all possible exchange terms, 

\F\) = \x)\exx...xx) + \x)\xex...xx) 

+ \x)\xxe...xx) + ... . (14) 

To measure the content of the beam, we introduce 
a "i/fi number" operator fi = This operator acts 

only on the beam states and gives (x\jl\x) = sin 2 a. To 
find the probability that the beam neutrino after inter- 
acting with N 2 background neutrinos will be measured 
as v^, we compute the expectation value of /2 in the final 
state of Eq. JT3J 

(F\ji\F) = a 2 (x\fi\x) 

x ((exx...x\ + (xex...x\ + ...) 

(\exx...x) + \xex...x) + ...) 
= a 2 sin 2 a((N% - N 2 ) cos 2 a + N 2 ) . (15) 

The last line is obtained by observing that in the 
sum there are N 2 "diagonal" terms of the type 
(x...xex...x\x...xex...x) = 1 and N 2 — N 2 "off-diagonal" 
terms of the type (x...xex...x\x...exx...x) = cos 2 a. 

The result in Eq. (|15f) shows how the connection be- 
tween a single scattering event and Eq. (JSJ) is made. For 
a single scattering event, N 2 = 1, one indeed finds that 



the conversion probability P Vb ^, v oc sin 2 a is maximized 
when the background consists of muon neutrinos. As the 
number of scattering events increases, however, the max- 
imum of conversion efficiency shifts to values of a < tt/2. 
In the limit of large 7V2, one finds P Uc _> u ^ oc N 2 sin 2 2a, 
precisely as predicted by Eq. (JSJ). 

We can now understand what happens to an electron 
neutrino propagating in a muon neutrino background. 
Eq. (|15f) shows that the conversion rate does not strictly 
vanish, as predicted by Eq. ijHJ), but is only proportional 
to N 2 , not to N 2 . As discussed before, the proportion- 
ality to N 2 is a feature of incoherent scattering. The co- 
herent scattering part is indeed strictly zero in this case: 
the states |e^^...) , \ fie fj.fi...), etc., are mutually orthog- 
onal, and hence have no overlap that could be coherently 
amplified. 

It can be readily seen that the part of the final state 
\F) that gets amplified contains the projection of the 
final background states onto the initial background state 
\xxx...). Decomposing the states in \F±) according to 

\x...xex...) — (x\e)\x...xxx...) + (y\e)\x...xyx...) , (16) 

where y is the state orthogonal to x, (y\x) — 0, we can 
write the result (|15|) as 

(F\fi\F) = a 2 (x\fi\x) 

x ((e\x)N 2 (xxx...x\ + 

(e\y)(yxx...x\ + (e\y)(xyx...x\ + ...) 

{(x\e)N 2 \xxx...x) + 

(y\e)\yxx...x) + (y\e)\xyx...x) + ...) 
= ^(^Al^XiVlK^le)! 2 + iV 2 |<2/|e>| 2 ). (17) 

We now return to the problem of two intersecting 
beams, one containing electron neutrinos and another 
containing neutrinos in the flavor superposition state v x . 
Suppose N\ neutrinos from the first beam interact with 
N 2 neutrinos from the second beam. As a result of the 
interaction, the system evolves over time St according to 

\eee...)\xxx...) — > \F) = \eee...)\xxx...) +ia\F\}, (18) 

where \Fi) contains a sum of terms with all possible ex- 
changes, 

\Fi) = (\xee...e) + \exe...e) + \eex...e) + ...) 

x (\exx...x) + \xex...x) + \xxe...x) + ...). (19) 

As in the derivation of Eq. (|17f) . to separate the coher- 
ent and incoherent parts, the states in \Fi) need to be 
projected on the corresponding initial states and orthog- 
onal states. This is done using Eq. (|16f) and 

\xee...) = (e|a;)|eee...) + ([i\x) \ fj.ee... ) . (20) 

As a result, we obtain 
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\Fi) = N 1 N 2 \(e\x)\ 2 \eee...)\xxx...) 

+ N 2 (/i\x}(x\e}(\[iee...) + \efj,e...) + ...)\xxx...) 
+ JVi(e|a;)(2/|e)|eee...)(|j/xa;...) + \xyx...) + ...) 

+ ( t i\x){y\e}(\fiee...} + \efte...) + ...)(\yxx...} + \xyx...) + ...). (21) 

I 



We observe that the term in the last line contains a sum 
of mutually orthogonal states. This term represents the 
part that will not amplify coherently, as can be seen by re- 
peating the arguments given in connection with Eq. 1)17(1. 
We further note that if one drops this term, at first order 

I 



Eqs. 1|22I23I24|) represent the central result of this sec- 
tion. They show that if we take the initial state to be 
a product of single-particle neutrino states and evolve it 
in time — carefully separating the coherent effects and 
dropping the incoherent ones — the final state obtained 
is again a product of single particle states, each rotated 
according to Eqs. ^23124(1 . No coherent superposition of 
many-particle states is formed. 

It is important to emphasize that to arrive at this con- 
clusion we only needed to consider elementary scatter- 
ing events (exchanges) between the beam and the back- 
ground. No assumptions, such as decoherence between 
mass eigenstates in the background, were necessary. 

I 



Eqs. (|23|l and Q24[l. describing the evolution of the one- 
particle neutrino states over small time St, then become 

I 



We observe that, although Eq. (|26|l was obtained for 
the initial states v e and v x , the derivation used only par- 



in a the final state equals 

\F) = \e'e'e'...e'e'}\x'x'x'...x'x'}, (22) 

where 



(23) 
(24) 

I 

IV. ONE-PARTICLE EVOLUTION EQUATIONS 

A. Evolution equation: active-active oscillations 



The construction in the preceding section can be used 
to obtain a differential equation describing the time evo- 
lution of a one-particle neutrino state. We first recall 
that in deriving Eq. I|21|l the flavor blind non-exchange 
interactions were omitted. It is easy to see that to in- 
clude their effect one should add an additional term, 
A^iiV2|eee...)|a;a;a;...), to the final state in Eq. (|2~T)> . 



(25) 



(26) 
(27) 



ticle exchanges and did not in any way rely on the par- 



|e') = \e)+zN 2 a[l/2x\(x\e)\ 2 \e) + (fi\x){x\e)\fx)}, 
\x') = \x)+iN 1 a[l/2x\(e\x)\ 2 \x) + {e\x)(y\e)\y)}. 



\Fi) = N 1 N 2 (l + \(e\x)\ 2 )\eee...)\xxx...) 

+ N 2 {^\x){x\e){\^ee...) + \efie...) + ...)\xxx...) 
+ JVi(e|a;)(2/|e)|eee...)(|ya;a;...) + \xyx...) + ...) 
+ {^\x){y\e){\^ee...) + |e//e...) + ...)(\yxx...) + \xyx...) + ...). 

I 



|e') = |e)+iJV 2 a[l/2x (1 + \(x\e)\ 2 )\e) + ^\x)(x\e)\^)}, 
\x') = |a;)+W 1 a[l/2x(l + |(e| a; )| 2 )| : r) + (e|a ; )(y|e)|y)]. 

I 
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ticular choice of the initial states. Therefore, at any time 
t, if the initial state of the beam neutrino is ip and the 
background neutrinos are all in the state 0, the evolution 
over a small time step is 



m + st)) - \m) = - 

[l/2x(l + |(^)| 2 )|V;) 



iV2G F N 2 (l - cos (3) /VSt x 
+ {^\4>){4>ty)\^)\, (28) 



where ij} is the state orthogonal to ip and we restored 
the coefficients, including the angular factor. 

Making use of the completeness relation \tp) + 
\4> ± )(iJj j -\ = 1, we rewrite Eq. as 



\1>(t + 6t))-W)) 
[1/2|V) + |^)(^)- 



= -iV2G F N 2 (l 
1/2|<«^>| 2 ^>]. 



cos 13) /VSt x 



(29) 



The evolution of a one-particle neutrino state is therefore 
described by the following differential equation: 



l : 

dt 


— H aa 


H aa 


= E 

i 

ri 



V2Gi 
V 



(1 -cos/? fo) ) x 



-\(^\^)\ 2 



(30) 



Here index (i) refers to a given particle for which the 
equation is written, and (j) runs over all other neutri- 
nos in the ensemble. The summation over the scattering 
angles (3^ was introduced to make this equation appli- 
cable to a more general case of neutrinos propagating in 
different directions. 



B. Evolution equation: active-sterile oscillations 

The preceding analysis dealt with flavor conversions 
between active neutrino states. The method developed 
there can be extended to describe the conversions be- 
tween active and sterile flavor states. In doing so, one 
should keep in mind two important differences between 
the two cases. In the active-sterile case: (i) only ac- 
tive components participate in the interactions, and (ii) 
the interaction amplitudes are proportional to the active 
content of the neutrino states involved. Performing cal- 
culations similar to those in Sees. IIIII and llV Al (see the 
Appendix) we find that the single-particle Hamiltonian 
for this case is given by 

B$ = £ -cos,**)) 

3 

x |(^)| e )| 2 [| e )( e |-l/2|^ W | e )| 2 ]. (31) 

As before, ip( % > denotes the state of the neutrino for which 
the evolution equation is written and cf>^' represents the 
flavor state of the jth neutrino in the background. It 



is worth noting that this Hamiltonian includes both the 
effects of the f-channel and the u-channel diagrams. 

For comparison, the standard Hamiltonian for an 
active-sterile neutrino system (the analogue of Eq. ©) 
is 

*<? = j: 2 ^(i-cos(3M)Hy j) i 

= E^p £ ( 1 - cos / 3fo) )K^' ) i e )i 2 i e )( e i-( 32 ) 



Here is the mixing angle of the jth neutrino in the 
background, cos 2 a^') = |(0^'|e)| 2 . 



C. Background of antineutrinos 

We now consider flavor transformation of the neutrino 
beam caused by the presence of antineutrinos in the 
background. For concreteness, let us envision a modi- 
fication of the thought experiment depicted in Fig. in 
which the beam contains neutrinos in a superposition 
of flavor states, v z = cos/3^ e + sin/3^, and the back- 
ground contains antineutrinos in a flavor superposition 
v x = cos aD e + sin a.v To fix the notation, let v w be the 
flavor state orthogonal to u z , (v w \v z ) = and v y be the 
flavor state orthogonal to v x , (D y \D x ) = 0. 

Consider what elementary processes are possible in this 
case. As can be easily seen, in addition to the t- and 
it-channel processes, v z v x — > v z v x and v z v x — > v x v Zl 
it is possible to have s-channel annihilation diagrams 
v z v z — ► v z v z ,v w v w ,v z v z ,T> w v w . Notice that only the v z 
component of the v x state participates in the s-channcl 
processes. 

The t-channel process is flavor-diagonal and hence does 
not cause flavor conversion, just like the corresponding 
process in the neutrino background. The w-channel pro- 
cess puts an antineutrino in the beam and, therefore, can- 
not interfere coherently with the incident neutrino wave. 
Any coherent flavor changes, therefore, can only be due 
to the s-channel process. 

Since only the v z component of the v x state partici- 
pates in the s-channel process, the beam neutrinos will 
not change flavor if the background antineutrinos are in 
the orthogonal flavor state, v w . While this is similar to 
what was found for a background of neutrinos in the state 
v w , on the microscopical level the two cases are quite dif- 
ferent. For the v w background, the amplitude of flavor 
conversion for a single elementary event is nonzero, but 
the conversion rate is only proportional to N, because the 
amplitudes add up incoherently. By comparison, for the 
v w background already the elementary amplitude van- 
ishes and the v w appearance rate in the beam is strictly 
zero. (Instead, D w antineutrinos will appear in the beam 
at the rate proportional to N, due to the u-channel pro- 
cess.) 
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The elementary scattering event can be written as 

\v z v x ) — ► \Vzp x ) - ia(v z \v x )(\v z v z ) + \v w v w )), (33) 

where the ^-channel process as well as the processes that 
put an antineutrino in the beam have been omitted. The 
minus sign appears because the amplitudes for neutrino- 
neutrino and neutrino-antineutrino scattering processes 
have opposite signs. 

As before, we project the final state on the initial states 
and orthogonal states, 

\v z v x ) — > \v z v x ) - ia{9 z \D x ) x 

{{v x \vz)\vgVx) + {v y \v z )\v z v y ) + 

v y )). (34) 

The rest of the argument proceeds in complete analoguey 
to the case of the neutrino background. The first three 
terms in parentheses in Eq. I)34|l will amplify coherently 
and the last term will not. Summing over many elemen- 
tary scattering events we obtain an expression similar to 
Eq. 1|21|) . which gives a one-particle evolution equation 



dt 



H a - a \lP 



cos/3 fe >) x 



l + \$ (J) )($ U) \-l\($ U) \1>®)\ i 



,(35) 



where <jyj) denotes the flavor state of the jth antineutrino 
in the background and the contributions of the i-channcl 
processes have been included. Thus, the effects of the 
neutrino and antineutrino backgrounds on the coherent 
neutrino flavor evolution have exactly the same form (but 
opposite signs!), even though at the microscopical level 
the two cases are different. 

When both neutrinos and antineutrinos are present in 
the background, their refractive effects add up linearly. 
Therefore, the Hamiltonian describing the flavor evolu- 
tion of a neutrino beam equals the sum of the two con- 
tributions, Eqs. (|3L)fl and (|35|l . Further generalization to 
include the effects of other matter (electrons, nucleons) 
and vacuum oscillations is straightforward. Just like in 
the case of the usual MSW effect, one should add the 
Hamiltonian terms H v&c and -ff ma t to the neutrino in- 
duced Hamiltonian (see Eq. QJ). 



D. Comparison to the standard results 

We now compare the one-particle Hamiltonians we 
have obtained, Eqs. (|30|l . I|31|) . and l|35|) . to the corre- 
sponding accepted results, Eqs. © and ipffil. We see 
that, while the accepted results are similar to ours, there 
are important differences: in all three cases, our Hamil- 
tonians contain additional terms. It is important to un- 
derstand both the origin of this difference and whether 
it leads to any physical consequences. 



First, we would like to establish whether the presence 
of the additional terms in our results affects the flavor 
evolution of the neutrino system. As can be readily seen, 
in all three cases the additional terms are proportional 
to the identity matrix in the flavor space. The evolution 
equations thus have the form 



^' = (#o + c(v,#r)v, 



(36) 



where Hq is the "standard" Hamiltonian given in Eq. © 
or 132fl and I is the identity matrix. We observe that, 
if ipo{t) is the solution of the equation itp' = Hoip, then 
-01 (t) given by 



ipi (t) = exp 



-i / c{Mt),Mt))dt 



Mt) (37) 



is a solution of Eq. (|36[) . This means that the extra terms 
give an overall shift to both energy levels, without affect- 
ing neutrino flavor evolution. 

Indeed, the physical origin of the difference in all three 
cases can be traced to the part of the interaction that 
does not change flavor. For concreteness, we for a mo- 
ment specialize to the first of the three cases, the active- 
active conversions in the neutrino background. It proves 
instructive to return to the evolution of the beam neu- 
trino over an infinitesimal time step. In our case, the 
result is given by Eq. (|26[l . while the accepted result, 
Eq. J2J), gives 

\e') = \e) + iN 2 a[(l + cos 2 a)\e) + sinacosa|/i)]. (38) 

The difference between the two is the factor of 1/2, which 
multiplies the state |e) (the flavor-preserving part) in the 
brackets in Eq. I|26|) . This factor, in turn, can be traced 
to Eq. I)25|l : it comes from the first term in |Fi), which 
must be split between the beam and the background to 
avoid overcounting. This is the origin of the factors of 
1/2 in Eqs. HO) and l(27j> . 

The situation is not unlike what happens in electro- 
statics. The interaction energy in a system of charges is 
given by 1/2 ^q,^ and the factor of 1/2 ensures that 
the interaction energy between pairs of charges is not 
counted twice. In our case, the extra terms serve the 
same purpose, to prevent counting the interaction energy 
twice. This can be seen as follows. Both active-active and 
active-sterile evolution equations are particular cases of 
a general case when the two flavor states have different 
couplings to the Z boson. As shown in the Appendix, the 
evolution equation in this general case can be written in 
a form JS7J: 



.d\il>« 
' dt 



H, 



5<* W 



(39) 



where is the generalization of the standard Hamil- 



tonian 
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(40) 



with G being the matrix of couplings, 



1 
T) 



(41) 



The second term in the evolution equation (|39|l has the 
form of the expectation value l/2(tp^\H^\ip^). This 
form makes explicit the physical meaning of this term 
as a correction to avoid double counting of the energy 
of the system. It is worth emphasizing that this term is 
always proportional to the identity in the flavor space, 
even when the two flavor states have different couplings. 
Physically, this is because the part of the interaction that 
needs to be split is always the part that conserves flavor 
(see the Appendix). 

The extra term, 1/2(iI>^\Hq\iP^), depends not only 
on the state of the background but also on the state of 
the beam itself. We therefore caution the reader that 
the superposition principle, which is always used for the 
MSW effect in normal matter, does not apply in this case. 
For instance, suppose that a beam neutrino which is v e at 
t = as a result of the evolution becomes a state v' and, 
similarly, a beam neutrino which is at t = becomes 
a state v" , Then, it is in general not true that the state 
v x = cosav e + sincw^ will become cosa^' + sinav" . 

It would be incorrect to conclude that the extra terms 
have no physical effect whatsoever. While they indeed 
do not change neutrino flavor evolution, they do mod- 
ify the absolute value of the refraction index of a neu- 
trino medium and, hence, at least in principle, change the 
bending of a neutrino beam in a dense neutrino medium 
with a density gradient. This effect is present even if 
there is only one neutrino flavor in the system. 



V. MORE ON THE ENTANGLED SYSTEM 

We have shown that if the neutrino system initially 
does not contain entangled states, such states are not 
formed as a result of coherent evolution in the system. It 
can be argued, moreover, that such evolution can lead to 
an effective loss of coherence between entangled states. 
To illustrate this, let us consider a beam in an entangled 
flavor state, 



|ent) = (| xxx. 



\yyy- 



•))/v / 2, 



(42) 



(here x and y are not necessarily orthogonal) propagating 
in the (uncntangled) background \zzz...). 

At time t = the expectation value of some opera- 
tor, for example, the number operator, ft, = J~l \/i)(n\, 



where the sum runs over all particles in the beam, is given 
by 

(cnt|/)|cnt) = -((xxx...|/t|xxx...) + (yyy...\ft\yyy...) 

+ (xxx...\ft\yyy...) + (yyy...\ft\xxx...)). (43) 

The first two terms on the right hand side simply count 
the muon neutrino content in the states |xxx...) and 
\yyy...) and the last two terms represent the effect of 
entanglement. 

Let us consider the effects of time evolution on the 
expectation value of ft. Each of the two terms in |ent) is 
a product of single-particle states and according to our 
earlier findings over time will remain a product of single- 
particle states. Let us write the state at time t = t\ as 



(\x'x'x'....)\z'z'z'...) + \y"y"y"....)\z"z"z"...))/V2, (44) 

where z', z", x', and y" are the results of solving a system 
of equations given in Eq. i|3U|) . (For example, z' and x' 
are found by solving the equations for the initial state 
\xxx...)\zzz...) .) 

The expectation value of ft in the state (|44|) is given 

by 

i((x'x'x'...|A|x'x'x'...) + {y"y"y"...W'y"y"-) 

+(x'x'x'...\ft\y"y"y"...)(z'z'z'...\z"z"z"...) 

+ (y"y"y"...\ft\x'x'x'...)(z"z"z"...\z'z'z'...)). (45) 

Since the states z' and z" will generically be different, the 
absolute value of the inner product (z"\z') will be < 1. 
The last two terms in Eq. I|45|l therefore contain sup- 
pression factors ((z"^')! 2 and vanish as the number of 
neutrinos in the background N2 is taken to infinity. As al- 
ready mentioned, these terms represent the entanglement 
between the two states; the system therefore behaves as 



if the beam was an incoherent mixture of jx'x'x'...) and 

\y"y"y"-)- 

Of course, rigorously speaking, the entanglement in- 
formation is not completely lost in the system. It may 
happen that at some time t the states z' and z" will 
be such that |(2r"|z')| = 1. In this case, the entangle- 
ment effect will reappear. We, however, regard this as 
an artificial arrangement and therefore maintain that for 
practical purposes the coherence is destroyed. 

It is curious to note that, as the entanglement effect 
reappears, the phases of the states \x'x'x'....)\z'z'z'...) 
and \y"y"y" ■■■■) \z" z" z" ...)) will have an effect on the ex- 
pectation value of ft and hence the phases due to the 
additional term, introduced in Eq. i|3U|) to avoid over- 
counting of energies, will have a physical effect. 
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VI. CONCLUSIONS 

In summary, we have developed a conceptually sim- 
ple, physical picture of coherent neutrino evolution in a 
medium of neutrinos and have shown how coherent inter- 
ference of many elementary scattering events gives rise to 
the refraction phenomenon. Unlike the case of ordinary 
matter, in which coherent scattering leaves the back- 
ground unchanged, in the neutrino background the scat- 
tering events change the state of the background, thus 
requiring a different type of analysis. We have found that 
only part of the elementary scattering amplitude is am- 
plified coherently This explains certain seemingly para- 
doxical results, such as why a neutrino traveling through 
the medium of neutrinos of opposite flavor does not un- 
dergo coherent flavor conversion. 

We have shown that refraction does not lead to the 
creation of entangled states in a neutrino system, i.e., if 
the state of the system is initially described by a prod- 
uct of single-particle states, the state remains a product 
of single-particle states as the system is evolved in time. 
Furthermore, the evolution effectively destroys initial en- 
tanglement in the system. It follows that for each neu- 
trino the result of the coherent evolution can be described 
by a one-particle Schrodinger equation, as was assumed 
a priori in the literature. 

We have derived the one-particle equation for active- 
active and active-sterile flavor transformation scenarios 
for a neutrino in a neutrino background. We also de- 
rived the equation for a neutrino in an antineutrino back- 
ground. In all these cases, we found that in order to avoid 
overcounting of the interaction energy one has to intro- 
duce an extra term in the evolution equations that is not 



present in the standard analyses. We have proven that 
this extra term does not affect the flavor evolution under 
normal conditions. It does, however, affect the value of 
the refraction index and hence the bending of a neutrino 
beam in a dense neutrino medium. 
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Appendix: the cases of active-sterile neutrinos and 
neutrinos with different couplings 

Let us derive the result (|31|l . which refers to a system 
of one active and one sterile neutrino, v e and v s . We 
start by considering the elementary interaction between 
two neutrinos, as discussed in Sec. Ill Bl The analogue of 
Eq. (J5J in this case is 



J 



dt 



I We{k)u e {p)) \ 
\v e {k)v s {p)) 
\v a (k)v e (p)) 
V \v s (k)v s (p)) J 



V2Gi 
V 



(1 - cos/3) 



/ 2 





V 



/ \ Ve {k)uM) \ 

\v e {k)v s (p)) 

Ws{k)Ve(p)) 

\ \v s {k)v s {p)) / 



(46) 



which shows that, given an initial state \v z (k)v x (p)) , only its active-active component, \v e (k)v e {p)) , is affected by the 
evolution: 



\v z {k)v x (p)) \v z {k)v x (0)) - idt^^IL(l - cos (3)(e\z)(e\x)\v e (k*)v e (p)) . 



(47) 



r 



Next, we apply this result to the case of two orthogonal 
neutrino beams, in the spirit of what was done in Sec. II I II 
We take neutrinos in the first and second beams to be in 
the states \z) and \x), respectively, and omit neutrino 
momenta for simplicity. Similarly to Eq. (|18f) . we get 



\zzz...) \xxx. 
where 



.) => \F) = \zzz...)\xxx...) +i2a\Fi), (48) 



\Fi) = (e\z) (e\x)(\ezz...z) + \zez...z) + \zze...z) + ...) 
x (\exx...x) + \xex...x) + \xxe...x) + ...) . (49) 



Notice that the effects of both the w-channel and t- 
channcl diagrams are included above. 



One can then follow the same procedure as in Sec. [H 
(see Eq. $2l\ ) and decompose |Fi) as follows: 
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\Fi) = (e\z}{e\x}[N 1 N 2 (z\e)(x\e}\zzz..}\xxx...} 

+ N2(w\e)(x\e)(\wzz...) + \zwz...) + ...)\xxx...) 

+ N 1 (z\e}(y\e}\zzz...)(\yxx...} + \xyx...) + ...) 

+ (w\e)(y\e)(\wzz...) + \zwz...) + ...){\yxx...) + \xyx...) + ...)] (50) 

where \y) and \w) are the orthogonal states to \x) and \z) respectively ((x\y) = (z\w) = 0). The last term in Eq. I|50|) 
is not coherently enhanced and therefore can be dropped. This allows us to obtain, at first order in a, a factorized 
form: 

\F) ~ \z'z'z'...)\x'x'x' ...) , (51) 

where 

\z') = \z) + iN 2 a\(e\z)\ 2 \{e\x)\ 2 \z) + i2N 2 a\(e\x)\ 2 (w\e)(e\z)\w) 

= |z)+ i 27V 2 a(e|^)|<e|x)| 2 |e)- l 7V 2 a|(e|z)| 2 |(e|x)| 2 |z) 
\x') = \x) +iJV 1 a|(e|z)| 2 |(e|a;)| 2 |a;) + i2N 1 a\(e\z)\ 2 (y\e){e\x)\y) . 

= \x) +i2N ia (e\x)\{e\z)\ 2 \e) -iN ia \{e\z)\ 2 \(e\x)\ 2 \x) . (52) 



From this a one-particle equation follows: 
.dj|0) =E 2^ (i _ cos/3W)) 

3 

x|(^)|e)| 2 [|e)( e |-l/2|(^|e)| 2 ] |^<*>) , 

(53) 

which proves the result l|31(l . 

We now discuss the generalization of our findings to 
two active neutrinos with different couplings to the Z bo- 
son. The study of this case provides a unified description 
of the results for active-active and active-sterile cases we 
have discussed so far. Furthermore, it allows us to com- 
pare our results with the corresponding discussion given 



inrcf. |lfj. 

Let us consider two neutrino states, v e and v p , and take 
v e as having the ordinary Standard Model coupling, g ei 
with the Z boson. The coupling of v p with the Z is defined 
as g p = rjg e . In terms of 77, the two-neutrinos equation, 



(|46p. is generalized by the replacements v s 



v, 



and 



2 









( 2 

7] 77 

7] 77 

\ 2t) 2 



(54) 



It is manifest that the active-active (Eq. ©) and active- 
sterile (Eq. H46f) ') cases are recovered for ij = 1 and 77 = 
respectively. The generalization of Eq. I|50|l has the form 



= N x N 2 {x\{z\B(ri)\z)\x)\zzz...)\xxx...) 

+ N 2 (x\(w\B(r])\z)\x}(\wzz...) + \zwz...) + ...)\xxx...) 

+ N 1 {y\{z\B{r 1 )\z)\x)\zzz...){\yxx...) + \xyx...) + ...) 

+ (y\(w\B(i])\z}\x}(\wzz...} + \zwz...) + ...)(\yxx...) + \xyx...) + ...) . (55) 
Once the last term in Eq. I|55(l is neglected, as discussed in Sec. IIIII one gets the effective one-particle equations: 

\z') = \z) + ^iN 2 a(x\(z\B( V )\z)\x)\z) + iN 2 a(x\(w\B{ V )\z)\x)\w) 

\x') = \x) + ^N ia (x\(z\B( V )\z)\x)\x) +iN 1 a(y\(z\B(r 1 )\z)\x)\y) . (56) 

I 



Notice that, as before, the factors of 1/2 arise as a re- Eq. (55} • 
suit of splitting of the first (flavor conserving) term in 
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Using the form Q54[l. and working out the matrix ele- where Hq has the form 
ments in Eq. I|56[l . it is possible to obtain the generaliza- 
tion of the result Ij53(l . We get 



.d\ip« 
1 dt 



H, 



|^«) , (57) 



H, 



(I) 



(58) 



Here the angular factors (1 — cos (3^) have been restored 
to account for general orientations of the neutrino mo- 
menta, and we use the matrix of the couplings (normal- 
ized to g e ) in the (y e , v p ) basis: 



1 

rj 



(59) 



Our result l|58(l coincides with the one-particle Hamilto 
nian given in ref. [Tol | . 

As an example, we give the explicit expression of H 
for the case of orthogonal beams with the neutrinos in the 



(») 



second beam being all in the same state, 



0') 



cos av„ 



sin av n 



V2G F N 2 

>< 

V 

2 „, i „,™2 , 



2 cos a + rj sin a 

rj cos a sin a rj cos 2 a + 2rj z sin z a 



rj cos a sin a 

.2 _ i r> 2 ;„2 



(60) 



It is easy to see that Eq. <|6(J|) reproduces the limiting 
cases of active-active {rj = 1) and active-sterile (77 = 0) 
neutrinos, Eqs. J7J) and (|3*2"1) . 
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